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Abstract 

We have extended, in most cases through 24th order, the series expansions of the dimer density 
in powers of the activity in the case of bipartite ((hyper)-simple-cubic and (hyper)-body-centered- 
cubic) lattices of dimensionalities 2 < d < 7. A numerical analysis of these data yields estimates of 
the exponents characterizing the Yang-Lee edge-singularities for lattice ferromagnetic spin-models 
as d varies between the lower and the upper critical dimensionalities. Our results are consistent 
with, but more extensive and sometimes more accurate than those obtained from the existing dimer 
series or from the estimates of related exponents for lattice animals, branched polymers and fluids. 
We mention also that it is possible to obtain estimates of the dimer constants from our series for 
the various lattices. 
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I. INTRODUCTION 



Let us consider a classical spin-1/2 Ising system with nearest-neighbor ferromagnetic 
interactions, in the presence of an external magnetic field H, on a finite <i-dimensional 
lattice of N sites. The partition function can be written as 

Z N (H, T) = J2 ex P [ K S S * S J + h X) s *] ■ C 1 ) 

conf <ij> i 

Here Sj = ±1 denotes an Ising spin variable associated to the site i of the lattice. The 
first sum extends to all configurations of the spins, the second to all distinct pairs < ij > 
of nearest-neighbor spins, the third to all spins. We have set K = J/kgT, with T the 
temperature, J > the exchange coupling of the spins, kg the Boltzmann constant, and 
h = H/kgT the reduced magnetic field. It is a classical result^ that, for real T, Z^{H,T) 
has a sequence of zeroes that can occur only on the imaginary axis of the complex h = 
h' + ih" plane. This statement remains true for a wide variety of models with ferromagnetic 
couplings, including Ising models of arbitrary spin^, n- vector models^"- (at least for n < 4), 
etc.. As long as the number N of sites remains finite, these zeroes give rise to a sequence 
of logarithmic branch-points of the free energy. For T above its critical value T c , when 
the thermodynamic limit N — > oo is taken, these cuts coalesce into two continuous lines of 
singularities along the imaginary axis of the /i-plane that extend to infinity and originate at 
the "Yang-Lee(YL) edges", h" = ±/i yi (T), with /i yi (T) > 0. The partition function must 
then be nonzero and the free energy must be analytic in a strip containing the real axis of the 
complex reduced field plane, so that no phase transition can occur for T > T c . As T — > T c , 
we have hy L (T) — > 0, so the gap between the endpoints of the singularity lines shrinks and 
they tend to pinch the real axis of the complex field plane giving rise, for T < T c , to the 
characteristic discontinuity associated with the spontaneous magnetization. 

At fixed T > T c , in the vicinity of a YL edge, the magnetization per spin is expected to 
exhibit a power-law behavior 7 

M(H,T)~\h±tti^ L (T)\° (2) 

controlled by an exponent a. 

The main features of these YL edge-singularities have been extensively studied. A priori 
the exponent a might depend on the temperature and the particular features of the model, 
however it was observed that it is universal, i.e. it depends only on the spatial dimensionality, 
but not on the structure of the lattice on which the spin model is defined. Moreover, 
due to the presence of the magnetic field, this exponent does not depend on the number 
of components of the spin. Hence for an n-vector ferromagnetic system, it must be the 
same as for an Ising system. Finally, the study of solvable models in one dimension or for 
d — > oo (the mean-field model) or, for any dimension, in the n — > oo limit (the spherical 
model), and finally the numerical study 7 of non-trivial Ising systems in d = 2 and d = 3 
dimensions, indicates that, for T > T c , the exponents a = a(d) are temperature- independent. 
Due to these universality features, the properties of the YL edges have been related^ to 
those of the ordinary critical points in the theory of second-order phase transitions. In 
this respect, an important step forward was to show that a renormalization-group approach 
to the calculation of the exponent a can be formulated^ in terms of an effective Landau- 
Ginzburg one-component scalar field theory with a cubic interaction and a purely imaginary 
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coupling. This entails the upper critical dimensionality d c = 6, at and above which the 
exponent takes its mean-field value a = 1/2. Moreover a can be expressed as an expansion 
in powers of e = 6 — d, whose coefficients have been so far computed^ through the third 
order. See Eq. (TT2~j) below. 

The interest of an accurate determination of a(d) is enhanced by the realization that the 
YL edge exponent is also related to other exponents characterizing the behavior of several 
quite different systems: 

a) The pressure for d- dimensional fluids, with repulsive-core interactions among their 
constituents, exhibits an unphysical singularity located at negative values of the activity 
and showing universal exponents^— <p(d) — cr(d) + 1. A similar property is true for lattice 
systems of hard squares (and higher- dimensional solids). 

b) The number-per-site of directed branched polymers with s bonds (or of directed site- or 
bond-animals^ 4 -^ of size s) on a <i-dimensional lattice, has a scaling behavior characterized 
by a power of s with universal exponent 1 ^— 0_d(gT) — cr{d — 1) + 1. The scaling behavior of 
the directed branched polymers has been mapped 1 ^ into that of the aforementioned fluids 
at the unphysical singularity, in one fewer space dimension. 

c) The universal exponents characterizing the behavior of the number-per-site of large 
isotropic branched polymers 2 ^ in a good solvent (or, equivalently 21 , of undirected site- or 
bond-animals) on a d-dimensional lattice can be expressed as (pi(d) = a(d—2)+2. The scaling 
behavior of the isotropically branched polymers has been mapped 2 ^ into that, mentioned 
above, of the fluids in two fewer dimensions. 

d) A field theoretic mode\M of the Anderson localization exhibits the same critical be- 
havior as the isotropic branched polymers mentioned above in c). 

Refs.fllli! provide us with tables of numerical estimates of these exponents for the 
systems of a), b) and c) in spatial dimensionalities < d < 8, which corroborate these expo- 
nent identifications. In particular, estimates of the exponent (f>i(d) for high- dimensional 
undirected site-animals have been obtained from 15th-order activity expansions 24 , and, 
more recently, have been complemented by extensive MonteCarlo simulations^ and exact 
enumerations 2 ^. 

We have recently calculated 2 ^ 1 ^, in most cases through 24th order, the high- 
temperature(HT) and low-field expansion of the magnetization in an external field for several 
models in the Ising universality class, including the conventional Ising spin-1/2 model, on a 
sequence of bipartite lattices of spatial dimensionalities 2 < d < 7. In this paper, we shall 
use these data, most not yet existing in the literature, to improve the accuracy and extend 
the range of direct estimates of the exponents of the YL edge-singularity. Our work is based 
on the observation-^— that this task can be accomplished restricting to the study of the 
Ising magnetization in a particular HT limit, since the exponents a(d) are temperature- 
independent. It was shown that, in this limit, for each value of d, the bivariate expansion 
of the magnetization in powers of the inverse temperature at low field reduces to a simpler 
univariate expansion of a quantity related to the density of a dimer system on the same 
lattice, in powers of a variable z that can be viewed as the activity of the dimers. Hence, 
the estimates of the exponents a(d) can be obtained from the simpler study of a singularity 
in the complex z-plane. 

The paper is organized as follows. In the second Section, we describe briefly the dimer 
model and sketch the arguments relating the Ising magnetization and the dimer density. In 
the third Section, we discuss the numerical estimates of the exponents a(d). In the fourth 
Section, we briefly mention the possibility of obtaining accurate estimates of the dimer 
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constants from the dimer series for several lattices. In the final Section, we summarize our 
work and draw some conclusions. 

TABLE I: Expansion coefficients d n of the dimer density per site p = XmLi d n z n for various 
lattices. In the case of the square lattice, the coefficients were already known 3 -^ through order 17, 
in the simple-cubic lattice case through order 15, in the body-centered-cubic case through order 
13. 
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II. DIMER DENSITY AND ISING MAGNETIZATION 

The hard dimer model in two space dimensions was introduced 3 ^ long ago to study the 
adsorption of diatomic molecules on a regular surface. On a d- dimensional lattice of TV 
sites it describes an assembly of linear objects each one covering a single bond of the lattice 
including the two endpoint sites, under the constraint of no partial or total overlap [hard 
dimers). The statistics of this model is described by the configurational macrocanonical 
partition function 

N/2 

~ n (z) = 1 + J29n(s)z s . (3) 

s=l 

Here Qn{s) denotes the number of ways of placing s dimers over the links of the iV-site 
lattice, and z is the dimer activity. 
In the thermodynamical limit 

oo 

E(z)= limlE^z)} 1 ^ = 1 + Y,9(s)z s (4) 

s=l 
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TABLE II: Expansion coefficients d n of the dimer density per site p = XmLi d n z n , in the case of 
the h4sc and of the h4bcc lattices, for which no data exist in the literature. 
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and the macrocanonical potential is denned as 

T(z)=M~(z)) = J2lsZ s . (5) 

s=l 

The dimer density per site p(z) is then expressed in terms of the potential T(z) as 

PW = ^ = E^ ( 6 ) 

s=l 

with c? n = ri7 n . Let us now summarize the simple argumenlr^r— relating the HT and 
low-field expansions of the spin-1/2 Ising model in the limit of infinite temperature with the 
low-density expansions for hard dimers on the same lattice. In the thermodynamical limit, 
the HT expansion of the Ising magnetization per spin can be written as 

00 1 1 

M(H, T) = c + 2(l-c 2 )J2jYl ^mic 2 ™- 1 (7) 

1=1 m=l 

with v = tanh(i^) and c = tanh(/i). In the elementary graphical representation of the HT 
and low-field expansion of the specific free energy, each term of order I in v and 2m in c is 
associated to a configuration of I bonds and 2m < 21 marked vertices on the lattice. The 
coefficients surviving in the particular limit in which v — > and c — > 00 with fixed z = vc 2 , 
are associated to configurations of I separated bonds whose endpoints are 21 marked vertices 



5 



TABLE III: Expansion coefficients d n of the dimer density per site p = YlnLi i n the case of 

the h5sc and of the h5bcc lattices, for which no data exist in the literature. 
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and therefore they are related to the number of ways of placing / hard dimers on the lattice 
in such a way that for the magnetization (i.e. the field-derivative of the free energy) one has 

M(H,T)/c^[l-2p(z)} (8) 

with 

00 

= (9) 
1=1 

the activity expansion of the dimer density. 

The behavior eq.([2]) of the magnetization nearby the YL edge must translate into 

p(*)~iz-zbr ( 10 ) 

for the dimer density, where 



z = - lim (tan(/i'' (T)) 2 tanh(X) (11) 

T— >oo 



is a real negative quantity. 



III. THE EXTENDED EXPANSIONS OF THE DIMER DENSITY 

Our recent calculation^!^, in most cases through order 24, of the HT expansion of the 
spin-1/2 Ising magnetization in the presence of an external field, for several bipartite lattices 
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TABLE IV: Expansion coefficients d n of the dimer density per site p = Y^=i d n z n , in the case of 
the h6sc and of the h6bcc lattices, for which no data exist in the literature. 
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of dimensionality 2 < d < 7, enables us to extend through the same order the expansions 
of the dimer density in powers of the activity z on these lattices. Such expansions were 
already tabulated in the literature^ in the case of the square(sq)(17 terms), the simple- 
cubic(sc)(15 terms), the body-centered-cubic(bcc)(13 terms) lattices. The known coefficients 
are reproduced by our calculation. Actually, a few more coefficients could be derived from the 
published HT and low-field expansions 3 -^ for the sq(19 terms) and the sc(17 terms) lattices. 
Also for the hyper-simple-cubic (h4sc) lattice in d = 4 spatial dimensions, 16 correct terms 
of the free-energy expansion were tabulated 35 in the same paper. However, these data were 
overlooked in Ref.[3l( and in the successive work. 

We have extended the activity expansions in the case of the h4sc and the hyper-bcc 
(h4bcc) lattices in d = 4 spatial dimensions through the 24th order. Our calculation for 
the higher- dimensional hyper-bcc lattices extends also through the 24th order, whereas in 
the case of the five- and six-dimensional hsc lattices, it reaches only the 22nd and the 21st, 
respectively. For dimension d = 7, our expansion on the hsc lattice extend only to the 20th 
order. These data were not yet available in the literature. 

The expansions for the h5sc, h6sc,...,h7sc lattices have been truncated at a slightly smaller 
order, because on the (hyper-)cubic lattices the computation of the embedding numbers for 
the highest-order graphs is too time-consuming to be completed in just a few days by a single 
personal computer. While in the case of the hbcc lattices the computational complexity 
of the calculation is independent^ of the lattice dimensionality, in the hsc case it grows 
exponentially with the space dimensionality. 

We have not yet extended the calculation of the series in the cases of non-bipartite lattices 
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TABLE V: Expansion coefficients d n of the dimer density per site p = Yl^=i d n z n , in the case of 
the seven-dimensional sc and bcc lattices, for which no data exist in the literature. 
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such as the triangular, face-centered-cubic and their high-dimensional generalizations. The 
coefficients of our extended expansions of the dimer density series for 2 < d < 7 are reported 
in the Tables H.. ED 



IV. ANALYSIS OF THE DIMER EXPANSIONS 

For the dimer density series, we have located the nearest negative singularity in the 
complex 2-plane and estimated its exponent a by the numerical tools currently used to 
extrapolate the behavior of power expansions to the border of their analyticity disk, in par- 
ticular employing the differential approximants (DAs). These are a generalization of the 
well known Pade approximant(PA) method that approximates the series under study by the 
solution of an ordinary inhomogeneous linear differential equation with polynomial coeffi- 
cients. For a detailed illustration of their definitions, of their properties, and for a general 
discussion of the uncertainties in the estimates of the singularity parameters determined by 
this method, we address the reader to Refs. 15.27.28 



Let us first consider the few known terms of the e-expansion of the exponent a 

/ ,n 1 1 79 o . 1 >/ , 10445 , „ n . 4x 

aid) = e e 2 + — ( 3 e 3 + O e 4 12 

v ; 2 12 3888 V 81 SV ; 1250712 ; v ; K J 

where ((n) denotes the Riemann zeta function. The expansion eq. f)12l) can be resummed 
following various procedures, all yielding consistent results for 1 < d < 6. For example, 
several simple continuous interpolations of a(d) with respect to the space dimensionality 
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d, can be founcr^i in the literature, which are obtained from low-order constrained Pade 
approximants. The continuous curve in FigJTJ is the plot vs d of a [3/2] Pade approximant 
formed from the e-expansion eq. (fl2"]) . that satisfies the additional constraints of reproducing 
the known exact values^ er(2) = —1/6 for d = 2 and <r(l) = —1/2 at the lower critical 
dimension d = 1. For d > 6, the exponent should retain the value a(d) = 1/2 and so this 
part of the curve is drawn flat. The constrained [3/2] PA has the expression 



The resummation by this approximant is visually indistinguishable from that obtained by 
the [6/0] PA subject to the constraint of reproducing; the exact results for d = 0, 1,2. The 
latter approximation was proposed in Ref.[12| to describe the behavior of the exponents 
of the unphysical singularity at negative activity for repulsive-core fluids as d varies. Al- 
ternative, but completely consistent interpolations can be obtained by other more refined 
resummation procedures^, including a Pade-Borel method and a conformal mapping tech- 
nique. For each dimension d, the small spread among the estimates of a(d) obtained by 
the various resummations methods of the expansion eq. (j!2p can be used to roughly limit an 
interval of possible values of a(d), whose width increases with e. We have reported these 
spreads in the eighth column of Table IVI1 

An e-expansion through the third order has been computed^ also for the leading 
correction-to-scaling exponent 9(d). Its values 9(2) = 0.83(1), 9(3) = 0.622(12), 0(4) = 
0.412(8) and 9(5) = 0.205(5) for the various dimensions have been estimated^ by con- 
strained PAs. Here we shall use these estimates without further discussion. 

From these studies of the e-expansion, we are thus led to expect that the exponent a(d) 
should lie between a = —1/2, its d — 1 value, and a = 1/2, the d = 6 value. For intermediate 
values of d, the exponent should take values not far from those suggested by the various 
resummation procedures. In conclusion, a(d) is expected to take small and possibly positive 
values for 2 < d < 6. The DAs are probably unable to measure accurately exponents in this 
range of values and therefore, as suggested in Ref.j3l[, it is more convenient to analyze the 
^-derivative of the dimer density, rather than p(z) itself, because it has a sharper singularity, 
characterized by the exponent a — 1. 

Our procedure of numerical analysis is the following. To begin with, we determine a 
central value zq of the location of the singularity by averaging over an appropriate sample of 
estimates from first- or second-order DAs, selected among those using most or all available 
series coefficients. Evident outliers are excluded from the sample. The uncertainty attached 
to zq is a small multiple of the spread of the reduced sample. Then we can evaluate also the 
corresponding exponent a(d) using a sample of second- or third-order DAs which are biased 
with the previously determined central value Zq, i.e. are built to be singular at Zq. This 
computational procedure allows only partially for the influence of the expected corrections 
to scaling on the estimates of the exponent, just because the DAs are more flexible than 
PAs. However the estimates so obtained still suggest a significant residual influence of 
the corrections to scaling, which is especially relevant in the case of the (hyper)-sc lattices 
and for large dimensionalities, and so we have tried to further improve our accuracy using 
also the well known proceduro^ 1 ^ of evaluating the exponent after performing the variable 
transformation 



on the series. By this simple prescription, the leading non-analytic correction to scaling of the 



a(d) 



0.00564924115e 3 - 0. 141867697c 2 + 0.390189052e + 1/2 



(13) 



-0.0852567479e 2 + 0.947044730e + 1 




(14) 
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TABLE VI: Estimates of the location zq and the exponent cr(d) of the nearest singularity of the 
dimer density in the complex activity plane z for several bipartite lattices of dimension 2 < d < 
7. In the first column, we have indicated the structure and dimensionality of the lattice. The 
second and fourth column show the results of Ref . jilT] . while the third and fifth contain the results 
of this paper. Notice that, for d > 5, our estimates for the hsc lattices are based on shorter 
series. In the sixth column, we have reported the values of cr(d) = cfi(d) — 1 obtained 12 from the 
exponent occurring for repulsive-core fluids. In the seventh column, we have shown estimates^ 5 - 
of the exponent cr(d) = <f>i(d + 2) — 2 derived from a high-precision MonteCarlo simulation of 
site-animals. In the eighth column, we have indicated ranges of estimates of a obtained from 
different resummation prescriptions of the e-expansion. The last column reports exact results for 
a, whenever known. The data in the last two columns depend only on the lattice dimensionality. 



lattice 


z (Ref. [31]) 


Zq (Our work) 


a (Ref. [31]) 


<T (Our work) 


Fluids (Ref. [12]) 


Animals (Ref. [25]) 


a e-exp.(Rcf.[9]) 


a Exact 


sq 


-0.088962(1) 


-0.0889642(2) 


-0.1645(20) 


-0.1662(5) 


-0.161(8) 


-0.165(6) 


0.146-0.166 


-1/6 


tri 


-0.056076(2) 




-0.1620(15) 












sc 


-0.052002(6) 


-0.0520268(2) 


0.096(6) 


0.077(2) 


0.0877(25) 


0.080(7) 


0.079-0.091 




bcc 


-0.037309(4) 


-0.0373125(2) 


0.090(5) 


0.076(2) 










fee 


-0.02413(18) 




0.097(11) 












h4sc 




-0.0365624(4) 




0.258(5) 


0.2648(15) 


0.261(12) 


0.262-0.266 




h4bcc 




-0.0170399(2) 




0.261(4) 










h5sc 




-0.0281849(4) 




0.401(9) 


0.402(5) 


0.40(2) 


0.399-0.400 




h5bcc 




-0.00814235(4) 




0.402(2) 










h6sc 




-0.0229462(6) 




0.460(50) 


0.465(35) 




1/2 


1/2 


h6bcc 




-0.0039832(2) 




0.475(30) 










h7sc 




-0.019360(3) 




0.495(8) 








1/2 


h7bcc 




-0.0019715(3) 




0.498(3) 











given series is approximately changed into an analytic one, which should be more accurately 
dealt with by the DAs. We observe that generally the estimates show little sensitivity to 
small variations of the bias value of Zq and of the correction exponents 6(d), so that our final 
uncertainties can generously allow also for the spreads of these input parameters. 

The estimates of the singularity location zq in the activity plane and of the exponent 
a(d) obtained by this kind of analysis are shown in the third and fifth columns of Table 
IVII and compared with those from a previous accurate DA analysis 31 of activity expansions 
shorter than those now available, that was performed along similar lines, but not explicitly 
allowing for the corrections to scaling. Our results generally lie within the range roughly 
suggested by the resummed e-expansion and are essentially consistent with the existing pure 
DA estimates^ 1 ^ for d < 4, but are more precise. It is worth to stress that our expansions 
of the dimer-density cover the whole range of dimensionalities between the lower and the 
upper critical dimension and therefore our results are also more extensive. For comparison, 
in the Table [VI] for each value of d, we have also reported estimates of a(d) obtained^ 
from the exponents (j)(d) associated to repulsive-core fluids and estimates obtained from the 
exponents 4>i(d+2) associated to the undirected site animals in d+2 dimensions. The former 
data were obtained^ by a DA analysis (which allowed for the corrections to scaling) of the 
14th-order fugacity-expansions for a model of a binary molecular mixture, with Gaussian 
interaction, in dimensions 1 < d < 6. The latter data were more recently obtained 25 by 
extensive MonteCarlo simulations performed in the case of (hyper)-simple-cubic lattices 
with 2 < d < 9. Both sets of estimates are essentially consistent with our results within 
their (sometimes larger) uncertainties. 
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It is also worth noticing that, if we allow explicitly for the corrections to scaling using the 
variable transformation defined by eq. (!14p . our estimates of <r(3) are somewhat smaller than 
the older results 31 from shorter dimer series, show a closer agreement with the constrained 
PA resummation of the e-expansion and exhibit more accurately the expected universality 
with respect to the lattice structure. The universality check is satisfactory also in d = 4. 
For d = 5 dimensions, the check is still as accurate, in spite of the small leading correction 
exponent, because of which other unaccounted higher corrections to scaling might also be 
significant. 

At the upper critical dimension d c — 6, it is expected that the power singularity of the 
dimer density is multiplicatively correctecU^* 3 -^^ by a power of a logarithm and that moreover 
also the leading corrections to scaling are logarithmic. Since the DAs are notoriously unfit 
to describe accurately this kind of singularity, it is not surprising that their convergence 
deteriorates. For lack of established and effective prescriptions to deal with this computation, 
we cannot do better than assigning to the estimate of 0"(6), an uncertainty much larger than 
the apparent spread of DAs. 

For completeness, we have computed a(d) also for the higher- dimensional hsc and hbcc 
lattices just to check that, as expected, above the upper critical dimension, the exponents 
a(d) retain the mean-field value, so that a(7) = 1/2, etc.. For graphical convenience, we 
have reported in the figure only the seven- dimensional result. 

In the Fig{TJ we have always reported the weighted average of the estimates for the sc- 
and the bcc-type lattices and attached to it a small multiple of the largest uncertainty of 
the single results. 



V. DIMER CONSTANTS 

Let us now show briefly that sufficiently long low-activity expansions of the dimer density 
can be of use also in estimating heuristically the constant hd that controls the exponential 
growth of the number of all possible dimer arrangements over the bonds of a A r -site d- 
dimensional simple-cubic lattice in the large N limit. This quantity is defined by 

h d = ln[ lim (5 Ar (^)) 1 / JV ] z=1 = T(l) (15) 

N— >oo 

and is often called monomer- dimer constant of the lattice under study, the term monomer 
referring to a site unoccupied by a dimer. 

It is more difficult to get from the series data good estimates also for the constant that 
characterizes the exponential growth of the number of the closely-packed dimer coverings 
(i.e. dimer arrangements such that every lattice site belongs to a dimer) for a iV-site d- 
dimensional simple-cubic lattice in the large N limit, 

h d = ]im[T(z) -p(z)\nz}. (16) 

2— >0O 

The quantity e hd is often called molecular freedom. 

The constants 41 hd and hd (which of course have nothing to do with the reduced magnetic 
field mentioned in the Introduction), are of interest in chemistry, combinatorial mathematics 
and information theory, but unfortunately their exact values are not known except^— in 
the case of I12. In particular, the quantities h,2, h%, and h§ are known rigorously only 
through their bounds. 
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Consider first the computation of hd- In the case of the sq lattice, a naive [12,12] PA of the 
^-expansion of V{z) yields the estimate h 2 ~ 0.662798.... This result agrees to six significant 
digits with the most precise determination 4 ^ 6 - h 2 = 0.6627989727(1) and is consistent with 
the tightest presently known bounds 0.66279897190 < h 2 < 0.662798972844913. Other 
accurate, but also non-rigorous estimates of h 2 can be found in the literature 2 ^ 4 ^ - — . In 
particular Ref.[50| reports a twelve figure determination of this constant. 

By the same quite simple prescription, in the case of the sc lattice, we find h 3 pa 0.7859..., 
a value which lies between the known 52 bounds: 0.7849602275 < h 3 < 0.7862023450. For 
the h4sc lattice, our PA estimate is ~ 0.880..., which should be compared with the lower 
bound 46 0.8638570485. These results show that very naive PAs yield rather accurate analytic 
continuations of the z-expansion of T(z) well outside its very small convergence disk. The 
precision of the estimates decreases as the space dimension d grows simply because the radius 
of convergence \z \ of the z-expansion of T(z) vanishes as d grows (see Table IVT1) . 

Unfortunately, this straightforward approach does not provide an analytic continuation 
of r(z) of acceptable accuracy for z » 1, that is necessary to compute also h d . We can try 
to mitigate this difficulty, performing in the z-expansion the simple change^ 1 ^ of variable 
z = t/(l — (q — l)t) 2 , where q is the coordination number of the lattice under study. As a first 
result, we can observe a drastic improvement in the apparent accuracy of the aforementioned 
estimates of the constants hd- For example, we obtain h 2 ~ 0.662798972.., which reproduces 
three additional digits of the best estimate of h 2 . For the three-dimensional sc lattice, the 
improved estimate is h% pa 0.7859660... and in four dimensions we obtain h 4 pa 0.88071788.... 
For the five- dimensional cubic lattice: h 5 pa 0.958123..., consistently with the known lower 
bound 0.94383303 < h 5 . Let us also stress that for the five- (and higher-dimensional) hsc 
lattices our expansions are slightly shorter. For example, they extend only through the 22nd 
order in the five- dimensional hsc case. 

Correspondingly, we have computed the values of the dimer density per site for z = 1 
obtaining: p 2 (l) ~ 0.31906155... (in complete agreement with the estimates of Ref. 45lJ47j ]). 
p 3 (l) pa 0.3421901..., p 4 (l) ~ 0.3579234... and p 5 (l) « 0.369580.... Other estimates for 
higher- dimensional hsc lattices and for the hbcc lattices can also be found in Table IVIH For 
d > 3, our estimates of h d for the hsc lattices compare well with the lower bounds provided 
by the inequality 4 ^ 1 ^ 

h d > -[-plnp-2(l -p)ln(l - p) + p\n2d - p] (17) 



with p = 4d+1 -f 5 + T . The meaning of the estimates of hd, summarized in Table IVHt is the 
following: we have reported only the decimal figures which appear to be stabilized in the 
three highest-order diagonal PAs that can be formed from the available dimer expansions. 

In terms of the new variable, the potential T(z) can be, with a better approximation, 
continued by PAs into a larger region of the complex z-plane, encompassing intermediate 
values of z such as z ~ 10. We have tried to infer the large z limit of the quantity T(z) — 
p(z)lnz in eq.f lTB"]) by observing that, in the intermediate z region, this quantity behaves 
linearly in 1/v^ to a fair approximation, and extrapolating this behavior to z = 00. We 
can thus conjecture estimates of the constants hd, that although generally reasonable and 
consistent with the known 46 bounds, unsurprisingly show an accuracy much lower than in 
the hd case and therefore will not be reported in the Table IVII1 The blame for this failure 
lays on the still insufficient extension of the region of accurate analytic continuation and 
the arbitrariness of the various extrapolation procedures, that make a safe evaluation of the 
uncertainties difficult. 
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TABLE VII: Our estimates of the dimer constants and the corresponding estimates of the dimer 
densities ^(l) (third column) for (hyper-)simple-cubic and (hyper-)body-centered-cubic lattices of 
various dimensionalities. In the second column we have reported lower bounds and in the fourth 
upper bounds for hd on the hsc lattice when available. No bounds are known for the hbcc lattices. 
In the text we have reported also estimates by other authors, when existing. 



hsc 


Lower limit 


This work 


Upper limit 


hbcc 


This work 


h 2 


0.66279897190 


0.662798972(1) 


0.662798972844913 






Pa(l) 




0.31906155(1) 








h 3 


0.7849602275 


0.7859660(1) 


0.7862023450 


^3 


0.8813479(1) 


Ps(l) 




0.3421901(1) 




P3(l) 


0.3584362(1) 




0.8638570485 


0.880718(1) 




Ha 


1.131880(1) 


Mi) 




0.3579234(1) 




Mi) 


0.392552(1) 


^5 


0.94383303 


0.95813(1) 






1.408080(1) 


Ml) 




0.369580(1) 




Mi) 


0.420248(1) 


he 


1.01129436 


1.023732(1) 




he 


1.703890(2) 


P6(l) 




0.37868106(1) 




Mi) 


0.4417731(1) 


h T 


1.06972606 


1.080759(2) 




h 7 


2.0142(2) 


Mi) 




0.38604998(2) 




Mi) 


0.457881(2) 



VI. CONCLUSIONS 

We have derived the coefficients of the series expansions of the dimer density in powers of 
the activity, in most cases through 24th-order, for several bipartite lattices of dimensionality 
2 < d < 7. Thus we have not only extended the existing series data, so far published only 
for d < 3, but also produced series for d > 4. An analysis of this set of data improves the 
accuracy in the determination of the locations and the exponents of the nearest real negative 
singularity of the dimer density in the complex activity plane. The exponents controlling the 
Yang-Lee edge-singularity of this class of ferromagnetic spin models in a wide range of space 
dimensions are related, in some cases by a simple dimensionality shift, with the exponents 
characterizing very different systems. The numerical estimates presently available for all 
these exponents show a good consistency among them and with the results of appropriate 
resummations of the renormalization-group e-expansion. Our dimer-density series are also 
shown to be of some use in estimating the dimer constants hd and hd- 
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